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Electric fields generate transverse flows near electrodes that sweep colloidal particles into densely packed 
assemblies. We interpret this behavior in terms of electrohydrodynamic motion stemming from distortions of 
the field by the particles that alter the body force distribution in the electrode charge polarization layer. A 
scaling analysis shows how the action of the applied electric field generates fluid motion that carries particles 
toward one another. The resulting fluid velocity is proportional to the square of the applied field and decreases 
inversely with frequency. Experimental measurements of the particle aggregation rate accord with the electro- 
hydrodynamic theory over a wide range of voltages and frequencies. 
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I. INTRODUCTION 

In its traditional form, the electrophoretic deposition of 
colloidal particles is analogous to sedimentation [1]. Since 
the electric field is a surrogate for gravity, (positively) 
charged particles move in the direction of the applied field. 
However, extensive experimentation during the last two de- 
cades has shown that particle motion transverse to the ap- 
plied electric field occurs near electrodes under a variety of 
conditions, as long as particles do not adhere to the surface 
[2-7]. As a consequence of the transverse motion, aggrega- 
tion takes place and the particles form densely packed planar 
clusters. Behavior of this sort is unexpected since the elec- 
trostatic interactions between identical particles — both Cou- 
lombic and induced dipole — are repulsive. Trau et al. [4,7] 
attributed the transverse motion to an electrohydrodynamic 
(EHD) mechanism wherein fluid flow carries particles to- 
wards one another. Particles near an electrode alter the local 
electric field and the action of these perturbations on the 
electrode polarization charge layer yields fluid motion. They 
showed experimentally that if the particles are monodisperse 
and spherical, the transverse movement can create two- and 
three-dimensional colloidal crystals in both dc and low fre- 
quency (cu< 10 kHz) ac fields. Bohmer’s observations of or- 
dering in dc fields [5] and the work of Yeh et al. [6] in dc and 
ac fields confirm the pervasive nature of the phenomenon. 
Aggregation of this sort has been observed with a variety of 
particles including nanoparticles [8,9] and biological cells 
[10,11], Nonspherical particles aggregate via EHD flow as 
well [9,10], though the aggregates may not always be crys- 
talline. Given its ubiquitous nature, the phenomenon should 
find diverse applications [12,13], 

Although it is generally agreed that electric-field induced 
fluid flow plays a central role in moving the particles to- 
gether, the detailed mechanisms remain somewhat unclear. 
Trau et al. ’s [4] EHD model is applicable in both ac and dc 
fields and consistent with experimental studies on homog- 
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enous electrodes as well as observations of particle motion 
toward ultraviolet-illuminated (high current) regions of semi- 
conductor electrodes [14]. Their model has also proved use- 
ful in understanding the formation of crystals with large in- 
terparticle separations [15-17], as well as superlattices in 
binary suspensions [18]. However, the theory was not elabo- 
rated to provide quantitative predictions of particle aggrega- 
tion behavior. 

To describe aggregation in dc fields, Solomentsev et al. 
[19] proposed an alternative mechanism based on electroos- 
mosis wherein flow due to the action of the steady field on 
the equilibrium diffuse layers around the particles draws 
them together. Since electroosmosis of this sort reverses in 
concert with the applied field and disappears at frequencies 
above a few hundred Hz, the model does not pertain to ag- 
gregation in ac fields. Sides [20,21] extended Trau et al .' s 
EHD model to investigate flow structure and the attendant 
particle movement in the vicinity of an electrode with a 
strong faradaic current, but concentration polarization at the 
electrode was ignored. 

In this work, we focus on aggregation of spherical par- 
ticles in fields oscillating at frequencies of several hundred 
Hz or higher. First, scaling expressions are derived to de- 
scribe the EHD flow engendered by field inhomogeneities 
near the polarization layer. The free charge density is mod- 
eled for a perfectly polarizable electrode with an alternating 
potential, neglecting the presence of particles. Then the point 
dipole approximation is used to describe perturbations in the 
field near the electrode due to the particles, taking account of 
mobile charge in the particle double layer. This produces 
lateral body forces in the polarization layer. Combining these 
ingredients yields a simple scaling expression for the EHD 
velocity that brings particles together. According to the 
analysis, the aggregation rate scales with the square of the 
field strength and inversely with frequency. We anticipate 
similar behaviors as to field strength and frequency with 
more complex particle shapes. 

To test the predictions of the scaling model, experimental 
techniques were developed to track the movements of a mul- 
titude of particles using video microscopy and image analy- 
sis. Then the disappearance of singlets (nonaggregated par- 
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FIG. 1 . A schematic diagram of the experimental apparatus (not 
to scale). A thin non-conductive spacer forms a cavity between two 
electrodes. The upper electrode is transparent — glass with a thin 
indium-tin-oxide layer — to facilitate optical microscopy. The lower 
silicon electrode is coated with a reflective layer of platinum. 


tides) was measured as a function of time over a wide range 
of field strengths and frequencies. The data were interpreted 
in terms of a second-order aggregation rate constant to es- 
tablish the effects of the main electric field parameters (ap- 
plied potential and frequency). Finally, comparison of the 
EHD model with experimental results shows that the scaling 
theory provides an accurate picture. 


II. SCALING THEORY 

To begin, consider the pure electrolyte between parallel 
electrodes (Fig. 1). When an electric potential is applied, 
counterions are attracted toward the electrodes and coions 
repelled. If electrochemical reactions are absent, the elec- 
trodes are “perfectly polarizable” and no current crosses a 
liquid-solid interface. As the applied potential varies in time, 
ion concentrations change to form and reform electrical 
double layers. Thus, charge on an electrode is balanced by 
charge in the diffuse layer. In the oscillating field, the fluid 
remains motionless and the pressure between the electrodes 
adjusts itself to maintain equilibrium. Solving the hydrostatic 
equation where pressure is balanced by electrical stresses 
with 


0=— VP + p if) E (1) 

and 

p ( ^ )= — ee 0 V 2 4> and E= — V cj> (2) 

shows that the pressure and the potential gradients are related 
as 


cillating at twice the imposed frequency. It follows as a cor- 
ollary that lateral disturbances to the uniform charge distri- 
bution cannot be balanced solely by pressure changes and 
flow ensues to provide the necessary viscous stress. 

The aggregation processes studied here arise from flow 
caused by the presence of small particles near an electrode 
but largely outside the charge layer. The particles perturb the 
otherwise uniform field, and the action of these perturbations 
on the charge in the polarization layer creates flow. It is 
stressed that the flow mechanism differs from classical elec- 
troosmosis in both the origin and dynamics of the charge 
distribution. In the classical model, the mobile equilibrium 
charge in the fluid balances charge that is chemically bound 
to solid-liquid interfaces. Here, the charge in the fluid fluc- 
tuates in response to the imposed, time-dependent potential. 
Because the induced charge is proportional to the applied 
field, electrical stresses scale on the square of the field and 
the response is justifiably designated as an electrohydrody- 
namic flow [22]. 

The relation between the tangential electrical stresses and 
velocity is revealed by inspecting the Helmholtz- 
Smoluchowski relation for steady electroosmosis in the equi- 
librium diffuse layer along a rigid, charged interface [23]. 
For the frequency range of interest, the charge layer thick- 
ness is characterized by the Debye scale, /c -1 [cf. Eq. (10)], 
and within the charge layer electrical stresses are balanced 
by viscous shear. For a 1-1 electrolyte, 

k~ 1 = \jsB 0 k B TI2n x e 2 . (4) 

The symbols stand for the bulk electrolyte concentration, 

; the charge on a proton, e\ Boltzmann’s constant, k B \ 
and the absolute temperature, T. In a 1-mM salt solution, 
k - 1 ~10 ran. Outside the Debye layer, the fluid velocity as- 
ymptotes to the Helmholtz-Smoluchowski value 


u 


* = 



( 5 ) 


Here, £ denotes the electrostatic potential at the electrode 
solid-fluid interface, /i is the fluid’s shear viscosity, and E, is 
the tangential component of the applied electric field. For a 
surface with £<0, the action of the field on the positively 
charged counterions produces fluid movement in the direc- 
tion of the field. It is instructive to rewrite this expression as 
a balance between the electric and viscous stresses on the 
Debye scale, i.e., 


p(z,t)=- l 2BE 0 (V(f>) 2 +C. (3) 

Here, the applied potential oscillates at a frequency w/2 ir, 
i.e., <p= <f> 0 (z)cos((ot); e is the fluid dielectric constant, and 
e 0 is the permittivity of free space. Since we are dealing with 
incompressible fluids, the constant of integration c can be 
taken as the pressure in the absence of the field, p 0 . Equa- 
tion (3) represents electrohydrostatic equilibrium and, since 
2 cos 2 («r)= 1 +cos(2o)t), the laterally homogeneous os- 
motic pressure arising from the action of the field on electric 
charge consists of a steady component plus a fluctuation os- 


u* ee 0 £ 

= — — E ,=qE,, ( 6 ) 

K K 

with Gauss’s law employed to relate the potential gradient 
and q, the total charge per unit area in the diffuse layer. 
According to Eq. (6), the induced velocity is proportional to 
the electrical stress per unit area, qE, . This proportionality 
will be used shortly to describe the velocity induced by per- 
turbations in the polarization charge. 

Further progress towards analytical expressions is facili- 
tated by linearization of the various conservation laws. First, 
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we calculate the charge in the polarization layer and then the 
EHD velocity. Hollingsworth and Saville [24] analyzed a 
linearized model of a dilute electrolyte between parallel 
plates with an applied potential difference, A 4>, oscillating at 
a frequency of u>I2tt. The potential is given by (the real part 
of) 


A <t> 


sinh( yz)(csch y) — iz yi' 2 coth y 
1 — i yv 2 coth y 


(7) 


where the various terms are defined as follows: 


A(f>=<f>\ z= - h -<f>\ z=h , 


y 2 = ( /c/i ) 2 



V 


2 _ 


(x) 

k 2 D 


Here, the electrode plate separation is 2h, with z = zlh , and 
D represents the ion diffusivity. 

For systems with Ka> 1 , particles of radius a maintain an 
average particle-to-surface gap that depends on a balance of 
several forces, including electrostatic repulsion, van der 
Waals attraction, hydrodynamic effects, and gravity. Our ob- 
servations show that particle movement in the electric field 
consists of Brownian motion superimposed on a directed 
translation. The rapidity of the Brownian fluctuations sug- 
gests that the particle-electrode gap is substantial. This is 
consistent with total internal reflection microscopy measure- 
ments by Fagan et al. [25], who found that the mean separa- 
tion distance is several Debye lengths, even in the presence 
of a strong oscillatory field. Outside the polarization layer, 
particles experience a uniform field directed perpendicular to 
the electrodes and from Eq. (7) this is 

_A (f> 
z=0 ~lh 

A</> 

( 8 ) 

Here, k denotes a unit vector in the z direction, cf. Fig. 2. 
Inasmuch as y 2 ~(*fi) 2 > 1 for w/2tt< 10 7 Hz, the field 
strength is mostly in-phase with the applied potential. 

To implement the scaling suggested by Eq. (6), we need 
the charge in the electrode polarization layer. In the fre- 
quency range of our experiments, 

/c 2 D~10 7 Hz>wl>/cD//i~400 Hz, (9) 

the free charge calculated from Eq. (7) is 



y csch y- i y v 2 coth y 
1 — i yv 2 coth y 


(a) 






(b) 



FIG. 2. Definition sketch depicting a spherical particle near an 
electrode, (a) The plan view depicting flow towards the test particle 
under conditions where the inequality [C' + (Dk 2 /o))C"]< 0 is sat- 
isfied. (b) An elevation view showing a particle located at x p outside 
the polarization layer; x t denotes a point in the polarization layer 
near the electrode, with x„-k=»0. 
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Of course the effects of suspended particles are also impor- 
tant but these are difficult to capture in detail without ana- 
lyzing the electrohydrodynamics of the region between an 
electrode and nearby particles. 

Next, to obtain the lateral component of the electric field 
due to a colloidal particle, we portray the particle’s influence 
in terms of a “point-dipole approximation.” The electric po- 
tential around a sphere in the uniform field well outside the 
diffuse layer (Fig. 2) is readily shown to be the real part of 



(x— x ) k— C„ 


(x-x p )-k 


(ID 


Here, E x = A<p/2h is the incident field strength, x and x p are 
vectors (scaled with the particle radius) for position and the 
location of the center of the particle, r is the (scaled) distance 
from the particle center, and C a is the (dimensionless, com- 
plex) “dipole coefficient.” Recall that for a dielectric particle 
in a dielectric medium, the dipole coefficient is a real quan- 
tity, given by the Clausius-Mossotti formula, viz., ( e p 
— e)/(e p + 2e) [26]. Here, we are interested in charged par- 
ticles suspended in an electrolyte where C„ is a complex 
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and the suspended particles will also influence matters, de- 
tailed numerical calculations will be useful to provide a com- 
plete quantitative picture of effects of ionic strength, particle 
size, etc. Nevertheless, the scale analysis furnishes an accu- 
rate picture of the effects of field strength and frequency. 

Equation (14) also provides a means to interpret the rate 
constant, k E , of a second order aggregation process in elec- 
trohydrodynamic terms (cf. Sec. III). If particles are carried 
towards one another by flow, then the flux towards a test 
particle across a circle of radius r is proportional to 
2nru(r)n(r), where n(r) represents the singlet particle 
concentration. Hence, the rate of aggregation per unit area 
will be second order in the number of particles per unit area 
with a rate constant whose magnitude behaves as 


FIG. 3. The dipole coefficient proportionality from Eq. (14) as a 
function of inverse frequency for a particle with a - 120-mV sur- 
face potential and 1-yum radius in 1-mM KC1. 


k E *=*2tra 


3ee 0 (A eft) 2 
flK 4h 2 


C' + 


Dk 2 
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(15) 


function of frequency and the properties of the particle 
(^-potential or surface charge and radius) and electrolyte 
(ionic strength, ion valences and mobilities, and dielectric 
constant). Accordingly, the dipole coefficient will be com- 
puted using the standard electrokinetic model [23], 

At x e , a point close to the electrode (near z=—h), the x 
component of the electric field around the dipole is the real 
part of 


E- i= 3E„C„ 


[(x t -x p )-k]x e -i _. a 


(Xp* k)(x e - i) 

■-3E„C 0 e 

r 


( 12 ) 


To establish the scaling of the EHD body force and its direc- 
tion, we proceed as follows. From Eqs. (10) and (12) and the 
properties of the trigonometric functions, the electric stress 
in the x direction at x c has a steady component proportional 
to 


qE,~ 3ee„ 


(A«/>) 2 
4 h 2 


Dk 2 

r' H C" 

O) 


(x p -k)(vi) 


(13) 


along with an oscillation at twice the imposed frequency. It 
follows that the tangential velocity due to EHD forces scales 
as 


3ee„ (A <f>) 2 
flK 4 h 2 


C' + 



(14) 


According to Eq. (14), the direction of the steady EHD mo- 
tion is established by the real and imaginary parts of 
C a (e £/ k b T , Ka ,a>a 2 / D , . . . ); flow will be towards the “test 
particle” if [C' + (Dk' 2 /co)C"]< 0. Since the dipole coeffi- 
cient is a weak function of frequency over this range, the 
velocity scales approximately as o> -1 . Figure 3 shows the 
outcome of a numerical calculation using the standard elec- 
trokinetic model [27]. Since the proximity of the electrode 


Thus, the rate of clustering is proportional to the square of 
the applied field {A<p/2h) 2 , and inversely proportional to the 
frequency a)/2n. Field strength and frequency are the main 
external control variables. As the experiments described in 
Sec. Ill demonstrate, this is the observed behavior. 

III. EXPERIMENTS 
Apparatus and materials 

The experimental cell consisted of two parallel electrodes 
separated by a 500- /u,m poly-dimethylsiloxane (PDMS) 
spacer (Fig. 1). The upper (grounded) electrode was glass 
with a thin coating of indium tin oxide, an optically trans- 
parent semiconductor (/? — 4—80, Delta Technologies); the 
lower silicon electrode was coated with a 150-nm-thick layer 
of platinum (Radiant Technologies). To clean the system, the 
electrodes and spacer were sonicated for 10 min in RBS-35 
detergent (Pierce Chemicals), thoroughly rinsed, sonicated 
again for 10 min in deionized water, and then dried with 
filtered nitrogen gas. Two batches of 2.1 -pm monodisperse 
polystyrene particles (Interfacial Dynamics Corp.) with dif- 
ferent surface charge densities were used. The particles were 
washed and diluted with the appropriate electrolyte solution 
to volume fractions of approximately 0.01. 

To begin an experiment, a suspension was placed in the 
cell and particles allowed to settle to form a dilute layer near 
the bottom electrode. Measurements of Brownian diffusion 
indicated that the particle-electrode separation was several 
hundred nanometers. Evidently the particles are kept out of 
contact with the electrode by a combination of electrostatic 
and steric effects. Oscillatory fields were then applied with a 
Tektronix function generator (PFG5505) and measured with 
a Tektronix digital oscilloscope (TDS2012). After acquiring 
a sequence of digital images at different stages in the clus- 
tering process, the voltage was turned off. Although Brown- 
ian motion breaks up particle clusters after aggregation it is a 
slow process; here a high frequency, high amplitude ac field 
was employed to expedite matters. Particle positions were 
observed with a reflection microscope (Mitutoyo) and re- 
corded using a CCD camera and a digital framegrabber. 
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Image analysis and kinetic interpretation 

Particle location was the main observable and, at low 
magnifications, 500-1000 particles could be imaged simul- 
taneously. We chose to analyze the behavior of a very large 
number of particles instead of focusing on a few binary 
events. Brownian motion often brings particles into close 
proximity and then separates them; concentrating on the be- 
havior of a few particles can be misleading. 

At the outset, the majority of the particles were isolated 
(i.e., not in apparent contact with other particles) and under- 
going Brownian motion. Upon application of the field, “sin- 
glets” moved together or toward larger aggregates to form 
clusters (Fig. 4). The observed motion was approximately 
two-dimensional and the rate of singlet disappearance served 
as a measure of the rate of aggregation induced by the elec- 
tric field. Since our focus is on the initial rate of particle 
aggregation, the crystallites are relatively small; it is clear 
from previous work that they will assemble into larger two- 
or three-dimensional crystals [4], 

To obtain the rate of singlet disappearance, images were 
recorded at one frame per second. Then standard image 
analysis techniques, including a thresholding step that sepa- 
rated individual particles (and higher order aggregates) from 
the background, were applied to each image. This yielded a 
list of clusters of various sizes and their locations. Although 
the use of low magnification complicates an accurate assess- 
ment of the number of particles in a given cluster, singlets 
have a relatively narrow distribution of “detected” sizes. A 
cluster that subtended p pixels (p min <p<p max ) was identi- 
fied as a singlet; clusters with p>p max were identified as 
higher order clusters. Inspection of a number of singlets in 
representative images indicated that the method is very ac- 
curate, with the majority of errors introduced at the edge of 
the image. Therefore, particles adjacent to the image edges 
were excluded from the analysis. 

Experiments at different field strengths indicate that the 
rate of singlet depletion is strongly affected by the character- 
istics of the electric field. To quantify the aggregation pro- 
cess, the rate of singlet disappearance can be expressed as 


dn ] 
dt 


= -2k xx n x n x 


N 


- 2 * 

i= 1 


\j*l*j* 


(16) 


where rij represents the concentration (per unit area) of clus- 
ters containing j particles, and k X j is a rate constant (with 
units of area/time) governing the rate at which singlets ag- 
gregate with clusters of size j. The factor of 2 accounts for 
the loss of two singlets each time a doublet forms; cluster 
breakup is neglected. To further simplify matters, we focused 
on the initial rate where particle disappearance was domi- 
nated by singlet-singlet events since higher order clusters 
were few in number. Under these conditions Eq. (16) simpli- 
fies to 



FIG. 4. Lateral aggregation in response to an electric field. Par- 
ticles are 2.1 -pm polystyrene in 1 .0-mM KC1. The applied field is 
500 Hz and 6 V, oriented out of the page, (a) Initially there are 407 
isolated particles (“singlets”) dispersed throughout the field of 
view, with a few additional particles in momentary contact due to 
Brownian motion, (b) After 10 s, the majority of particles have 
aggregated into clusters, with 87 singlets remaining, (c) After 30 s, 
only 19 singlets remain. 


— = 1 + k E n° x t. (18) 


SO 


dn '__. 2 
dt - k£n >’ 


The factor of 2 has been incorporated into the rate constant, 
k E , and n x is the initial particle concentration. According to 
Eq. (18), the inverse singlet concentration increases linearly 
with time, with slope equal to k E n x . 
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FIG. 5. Singlet concentrations (normalized on the initial concen- 
tration) as a function of time for different applied voltages at a 
frequency of 1 kHz. Squares, 6 V; downward triangles, 5 V; circles, 
4 V; diamonds, 3 V; upward triangles, 2 V. Upon application of the 
field, the singlet concentration decreases as aggregation ensues. Par- 
ticles are 2.7-/xm polystyrene in 1.0-mM KC1. The initial singlet 
concentration for each voltage is approximately one particle per 
70 /j, m 2 (~ 14 000 particles/mm 2 ). 

Aggregation kinetics 

By carrying out the analysis procedure on a sequence of 
images collected during an experiment, the number of sin- 
glets in a fixed viewing area was determined as a function of 
time. Representative data are shown in Fig. 5. To ensure that 
singlet-singlet events predominated, data analysis terminated 
when 20% of the singlets remained or 30 s after the field was 
applied, depending on which occurred first. As the data show, 
before application of the field, n x remains relatively constant, 
with small fluctuations as particles momentarily encounter 
one another due to Brownian motion. Upon application of 
the field, n , decreases rapidly with a rate dependent on the 
applied potential. Once the field is removed, the number of 
singlets increases as the aggregates disperse by Brownian 
motion (data not shown). Rescaling the data as «"/«! (Fig. 6) 



'<*) 


FIG. 6. Scaled (inverse) singlet concentration data from Fig. 5. 
Symbols are the same as in Fig. 5. Solid lines are linear fits to the 
data; the slopes of the lines yield the aggregation rate constants (see 
the text). 



FIG. 7. Aggregation rate constants obtained from a series of 52 
sequential experiments, each at 500 Hz and 4 V. A histogram (inset) 
shows the number of experiments with a given rate constant. 

shows that the inverse concentration increases linearly with 
time, consistent with Eq. (18). 

To test reproducibility, a sequence of 52 experiments was 
carried out with a 500-Hz, 4-V applied potential. As Fig. 7 
indicates, the rate constants are quite reproducible. A histo- 
gram (Fig. 7, inset) suggests a quasi-Gaussian distribution 
with a standard deviation of (approximately) 10% of the 
mean. Preliminary experiments indicated that one cause for 
the scatter is variation in the initial particle configuration. If, 
on average, particles are much closer to one another (i.e., 
higher concentrations) at the onset of the electric field, then 
the apparent rate of aggregation will be correspondingly 
larger due to higher order encounters. To minimize this ef- 
fect, the particles were moved into a widely separated con- 
figuration prior to each experiment using a 10-kHz, 10-V 
applied field. This provided a reproducible random initial 
condition for each experiment. However, it was difficult to 
precisely control the initial singlet concentration, which var- 
ied between experiments due to random aggregation and de- 
aggregation via Brownian motion. This partially accounts for 
the scatter in the measured rate constants, since contributions 
from higher order clusters are neglected in the analysis. 

To confirm that the particles were randomly positioned at 
the start of an experiment, the degree of disorder was quan- 
tified using a well-defined “translational order measure,” T\ 
T= 0 for perfect disorder and T = 1 for perfect order [28]. In 
the set of 52 experiments, the average value of T was 
0.024± 0.014, indicating that the ac field procedure was ef- 
fective in arranging the particles into disordered initial con- 
figurations. 

IV. COMPARISON OF THEORY AND EXPERIMENTS 

To probe the validity of the scaling expression set out in 
Eq. (15), the aggregation rate was measured as a function of 
the field strength and frequency. Rate constants (extracted 
from the data illustrated in Fig. 6) are plotted against the 
square of the applied voltage amplitude for frequencies of 
500 and 1000 Hz (Fig. 8). Ten experiments were performed 
at each amplitude and frequency. The aggregation rate con- 
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FIG. 8. Aggregation rate constant-applied potential relation. 
Particles are 2.7-yum polystyrene in 1.0-mM KC1. Each point is the 
mean of ten separate experiments; error bars represent one standard 
deviation. Two sets of experiments were done at 500 Hz with par- 
ticles of same size but different electrophoretic mobilities: circles, 
M-— 8.6 yumcm/V s; squares, M=— 9.5 yumcm/V s. Another set 
of experiments was carried out at 1000 Hz (triangles). Solid lines 
were derived by linear regression. 

stant increases linearly with (A <f>) 2 , in agreement with the 
scale analysis, and this dependence extends over an order of 
magnitude at both frequencies. The (A cf>) 2 dependence was 
replicated with a different particle, albeit with a similar £ 
potential (Fig. 8). It is reassuring to note that the numerical 
magnitude of the measured rate constant k E *= 10~" m 2 /s 
(Fig. 8), is consistent with that estimated using Eq. (15). 

Another series of experiments was employed to investi- 
gate the effect of frequency (Fig. 9). For polystyrene par- 
ticles in 1.0-mM KC1, the rate constant changed as 1/to for 
frequencies larger than 1000 Hz, in agreement with Eq. (15). 
The solid lines are correlations for frequencies greater than 1 
kHz, chosen because it is an order of magnitude larger than 
the minimum frequency kD/H~ 100 Hz for which the scal- 
ing analysis is valid. The inverse frequency dependence 
found here is similar to that identified by Kim et al. [29], 
who found a to -12 dependence between 30 and 500 Hz for 
10-/u,m polystyrene spheres in 0.1-mM NaHC0 3 . A key dif- 
ference is that we observe a maximum aggregation rate at 
— 500 Hz, with the rate constant decreasing as <o— >0 (data 
not shown). Since the inverse frequency dependence is only 
valid for frequencies much greater than — 100 Hz, it is not 
surprising that slightly different behavior is observed at 
lower frequencies. One explanation for diminished aggrega- 
tion rates at low frequencies may be that significant vertical 
electrophoretic motion is present; such motion was clearly 
observable as particles moved slightly in and out of the focal 
plane. Electrochemical reactions may also play a role at 
lower frequencies [20,25]. At higher frequencies, induced 
dipole-dipole repulsion tends to inhibit aggregation [15-18], 
and presumably affects the maximum frequency at which 
particles aggregate (i.e., the k E = 0 intercepts in Fig. 9). 


FIG. 9. Aggregation rate constant-frequency relation for two 
applied potentials. Symbols: circles, A(/> = 4 V; squares, A</> 
= 6 V. Each data point is an average from ten experiments; error 
bars are one standard deviation. Solid lines are linear fits for fre- 
quencies of 1 kHz and larger (inverse frequencies of 1 kHz -1 or 
smaller). 

Nonetheless, over the range of frequencies where electro- 
phoretic motion is negligible and the free charge is given by 
Eq. (10), the rate of aggregation decays inversely with fre- 
quency. 

V. CONCLUSIONS 

Scaling expressions for the EHD fluid velocity engen- 
dered by the dipole field of a polarized particle were derived 
and it was shown that the direction of flow depends on the 
sign of the dipole coefficient. According to the analysis, the 
velocity is proportional to the square of the applied potential 
and inversely proportional to the frequency. To test these 
scaling predictions, the behavior of large numbers of par- 
ticles was tracked and the aggregation kinetics quantified in 
terms of a second-order rate constant. The experimental re- 
sults are all consistent with the scaling analysis. 

Although the analysis presented here focused on oscilla- 
tory fields and monodisperse suspensions, the approach is 
applicable in other situations. For example, EHD flow occurs 
in steady fields, although analysis is complicated by several 
additional effects (e.g., electrochemical reactions and elec- 
troosmotic flow [19]). In addition, the scale analysis can be 
extended to binary or ternary suspensions [18] although 
more detailed models will be required to describe such sys- 
tems. 
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